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The KEK to Kamioka long-baseline neutrino experiment (K2K) is the first accelerator-based 
experiment with hundreds of km neutrino path length. K2K focuses on the study of the 
existence of neutrino oscillations in disappearance that is observed in atmospheric neutrinos. 
With nearly pure muon neutrino beam produced with 5.6 x 10 19 protons from 12 GeV PS 
on target, 56 neutrino events have been observed at Super-Kamiokande, the far detector at 
250 km distance. The expectation is 80.6lgQ derived from the measurement at the near site. 
The experiment expects to accumulate 10 20 protons on target, providing sufficient statistics 
to study neutrino oscillations by spectral analysis for z/ M disappearance. 



Introduction 



The existence of neutrino oscillations implies that neutrinos are massive and that flavor and 
mass eigenstates are mixed in the lepton sector, too. Experiments on atmospheric neutrinos 
have found a significant deficit in the flux of Vu which have travelled an earth scale distanc(0'§i. 
The interpretation of these results provides not only strong evidence of — ► v T (or — > v s tertte) 
oscillationiaO but also evidence for different mixing properties in the lepton and quark sectors. 
For a neutrino energy E^GeV) and a distance from the source L(km), the oscillation probability 
can be written in terms of the mixing angle 9 and the difference of the squared masses Am 2 



(eV ) in two flavor approximation as P(vu 



sin 29 sin 2 (1.27 Am 2 L/E v 



The KEK to Kamioka long-baseline neutrino experiment (K2K) is the first accelerator- 
based experiment with hundreds of km neutrino path lengtra. The intense, nearly pure neutrino 
beam (98.2% !/„, 1.3% v e , and 0.5% V a ) has an average L/E v w 200(L = 250km, (E„) w 



1.3GeV). The neutrino beam properties are measured just after production, and the kinematics 
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of parent pions are measured in situ to extrapolate the measurements at the near detector to the 
expectation at the far detector. K2K focuses on the study of the existence of neutrino oscillations 
in Vn disappearance that is observed in atmospheric neutrinos, and on the search for z/„ — ► v e 
oscillation with well understood flux and neutrino composition in the Am 2 > 2 x 10 _3 eV 2 
region. In this paper, the results on the event rates with data taken from June 1999 to July 
2001 corresponding to 5.6 x 10 19 protons on target are described. 



2 Experimental Setup and Analysis 

The K2K neutrino beanfli is a horn-focused wide band beam. Figure [l] shows the setup 
at KEK site. The primary beam for K2K is 12 GeV kinetic energy protons from the KEK 
proton-synchrotron (KEK-PSp. Every 2.2 s, approximately 6 x 10 12 protons in nine bunches 
are fast-extracted in a single turn, making a 1.1 fis beam spill. These protons are focused onto 
a 30 mm diameter, 66 cm long aluminum targetwhich is a current carrying element in the first 
of a pair of horn magnets operating at 250 kMQ. This design maximizes the efficiency of the 
toroidal magnetic field to focus positive pions produced in proton-aluminum interactions, while 
sweeping out negative secondary particles. 




Figure 1: Near site setup at KEK. 



Downstream of the horn system, before the 200 m long decay volume where the pions decay 
to Ufj, and muons, a gas-Cherenkov detector (PIMONjU is occasionally put in the beam to 
measure the kinematics of the pions after their production and subsequent focusing. After the 
decay volume, an iron and concrete beam dump stops essentially all charged particles except 
muons withan energy greater than 5.5 GeV. Downstream of the dump there is a "muon monitor" 
(MUMONp consisting of a segmented ionization chamber and an array of silicon pad detectors. 
This monitors the residual muons spill-by-spill to check beam centering and muon yield. The 
ionization chamber consists of 5 cm wide strips covering roughly 2 x 2m 2 area with separate 
planes for horizontal and vertical read-out. The 26 silicon pads, each of ~ 10cm 2 transverse 
area, are distributed through a 3 x 3m 2 area. 

Neutrino interactions near the production site are measured by a set of detectors with com- 
plementary abilities. The detectors are located approximately 300 m from the pion production 
target with approximately 70 m of earth eliminating virtually all prompt beam products other 
than neutrinos. A one kiloton water Cherenkov detector (lkt) uses the same technology and 



analysis algorithms as the far detector, Super-Kamiokande (SK). It has 680 20" photomulti- 
plier tubes (PMTs) on a 70 cm grid lining a 8.6m diameter, 8.6m high cylinder. The PMTs 
themselves and their arrangement are the same as in SK, giving the same fractional coverage 
by photo-cathode (40%). A scintillating fiber detector (SciFiP with a 6 ton water target has 
tracking capability, and allows discrimination between different types of interactions such as 
quasi-elastic or inelastic. Downstream of the SciFithere is a lead glass array for tagging electro- 
magnetic showers. A muon range detector (MRD^J, measures the energy, angle, and production 
point of muons from charged current (CC) interactions. It covers 7.6 x 7.6m 2 transverse area 
with four 10-cm iron plates, followed by eight 20-cm iron plates, all interleaved with drift tubes. 
The total mass is 915 tons. The 6632 drift tubes each has 5x7cm 2 cross section, and are arranged 
in horizontal and vertical read-out planes. 

The far detector for K2K is the SK detector located in the Kamioka Observatory, Insti- 
tute for Cosmic Ray Research (ICRR), University of Tokyo, which has been taking data since 
1996. Its performance and results have been documented in the literature! Event selection for 
this detector uses timing synchronization with the KEK-PS via the Global Positioning System 
(GPS0. 

The beam-line was aligned by GPS position surveji. The precision of this survey for the 
line from the target to the far detector is better than 0.01 mr and the construction precision for 
the near site alignment is better than 0.1 mr. The predicted neutrino spectrum at 250 km is 
approximately the same over nearly 1 km, giving a required accuracy of 3mr for the pointing of 
the beam. 

The steering and the monitoring of the v„ beam are carried out based on the MUMON 
measurement. The ia, and muons in the decay volume originate from the same pion decays, 
so the measured center position of the muon profile is correlated with the v„ beam direction. 
The r.m.s. muon profile width is ~ 90cm for both the horizontal and vertical directions. These 
measurements are used in tuning the beam direction at the start of every beam period, usually 
once per month. During beam tuning, the primary proton position at the target is adjusted so 
that the muon profile is centered on the SK direction within 0.1 mr. The center of the muon 
profile is measured to be stable spill- by-spill within ±lmr in each dimension throughout the 
whole experimental period. The stability of the muon yield is directly related to the stability 
of the horn-focused za, beam. The yield normalized to proton intensity measured by a current 
transformer is stable spill-by-spill within a measurement uncertainty of 2%. 

The characteristics and stability of the z^ beam itself are directly monitored at the near 
site by the MRD using interactions with iron. The large transverse area of MRD makes 
it possible to measure the beam direction and width (profile), and the large mass makes it 
possible to measure the time stability of the za, event rate, profile and spectrum. The location 
of the center of the profile gives the beam direction. For the profile measurement, the starting 
point of a reconstructed muon track is regarded as the vertex of the z^-iron interaction. The 
vertex distribution is corrected for geometrical acceptance to obtain the beam profile. The 
direction has been stable within ±1 mr throughout all data-taking periods. The measured rate 
of z/^-iron events is typically 0.05 events/spill and can be monitored on a daily basis with good 
statistics. The event rate normalized to proton intensity is stable within the statistical error 
of the MRD measurement. The muon energy and angular distributions are also continuously 
monitored. They show no change as a function of time, implying the energy spectrum is 
stable throughout this period. 

Comparison between expected and observed numbers of z^ events at the far site is done with 
this knowledge of measured beam stability. To predict the z/„ beam characteristics at the far site, 
a normalization measurement at the near site and the extrapolation of the information from near 
to far are necessary. For the rate normalization, the lkt is used so that any detector or analysis 
bias is suppressed. For the extrapolation, the beam simulation is used, which is validated by the 



PIMON measurement as described below. This simulation is based on GEANT0 with detailed 
description of materials and magnetic fields in the target region and decay volume. It uses as 
input a measurement of the primary beam profile at the target. Primary proton interactions on 
aluminum are modeled with a parameterization of hadron production dat£l. Other hadronic 
interactions are treated by GEANT-CALOF0. 

Once the kinematic distribution of pions after production and focusing is known, it is possible 
to predict the i/„ spectrum at any distance from the source. The pion momentum and angular 
distribution is measured by the PIMON. Cherenkov photons generated by pions in the detector 
are focused by a pie-shaped spherical mirror to an array of PMTs of 8 mm effective diameter 
in its focal plane. Photons emitted by particles with the same velocity and incident angle with 
respect to the mirror arrive at the same position in the focal plane. This is independent of 
photon incident position with respect to the mirror. The photon distribution on the PMT 
array is a superposition of slices of the Cherenkov rings from particles of various velocities 
and angles. Measurements are made at seven indices of refraction, controlled by gas pressure, 
to give additional information on the velocity distribution of the particles. The indices of 
refraction are chosen so that the corresponding pion momentum interval is 400 MeV/c. The 
pion two-dimensional distribution of momentum versus angle is derived by unfolding the photon 
distribution data with various indices of refraction. To avoid background from 12 GeV primary 
protons, the index of refraction is adjusted below the Cherenkov threshold of 12 GeV protons 
(1 — (3 ~2.6 x 10 -3 ). As a consequence, analysis is done for P n > 2GeV/c, giving the energy 
spectral shape above 1 GeV. Figure || shows the inferred energy spectral shape at the near site 
and the far to near flux ratio along with the beam simulation result. The beam simulation 
is well validated by the PIMON measurement without any tuning. 
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Figure 2: The top figure shows the energy spectral shape at the near site and the bottom figure shows the 
far to near flux ratio. The histograms are from the beam simulation results. The data points are derived from 
the PIMON measurement. For the top figure, normalization is done by area above 1 GeV. The vertical error bars 

for the data points reflect the total uncertainty. 



The interaction rate at the near site is measured in the lkt by detecting Cherenkov 
light emitted from produced charged particles. The vertices and directions of Cherenkov rings 
are reconstructed with the same methods as in In the lkt, however, an analog sum of 

signals from all PMTs is recorded by an FADC to select beam spills with only one event. The 



definition of events per spill is the number of peaks of the FADC signal with > 1000 collected 
photo-electrons (p.e.) (~100MeV deposited energy). The average number of events observed in 
the full detector is about 0.2/spill including entering background events due to cosmic rays or 

induced muons from upstream (about 1.5% of events in the fiducial volume defined below). 
Thus > 2 events occur in about 10% of those spills with > 1 event. Neutrino interactions are 
selected by requiring: a) There is no detector activity in the 1.2/xs preceding the beam spill, b) 
Only a single event is observed in the FADC peak search for that spill, c) The reconstructed 
vertex is inside the 25 t fiducial volume defined by a 2 m radius, 2 m long cylinder along the 
beam axis, in the upstream side of the detector. The detection efficiency of the lkt detector 
for detecting CC interactions is 87% and for neutral current inelastic (NCi ne i) interactions it is 
55%. The expected ratio of CC to NCi ne i gives an overall efficiency of 72%. The main source 
of inefficiency is the 1000 p.e. threshold of the peak search in the FADC signal. The neutrino 
interaction simulator used for efficiency calculations is the same as that used in all SK analyses. 
Finally the number of events in lkt is corrected for spills with multiple events. The average 
v„ event rate per proton hitting the target is 3.2 x 10 -15 . Correcting for efficiencies, relative 
target masses, and detector live times, the expected signal at the far site is estimated by 
applying the extrapolation from the experimentally validated beam simulation. 

Accelerator-produced neutrino interactions at the far site, SK, are selected by comparing two 
Universal Time Coordinated (UTC) time stamps from the GPS system, Tkek for the KEK-PS 
beam spill start time and T^k for the SK trigger time. The time difference between two UTC 
time stamps, A(T) = T^k — Tkek ~ TOF where TOF is the time of flight of a neutrino, should 
be distributed around the interval from to l.lyus to match the width of the beam spill of the 
KEK-PS. Since the measured uncertainty of the synchronization accuracy for the two sites is 
< 200ns, beam-induced events are selected in a 1.5//S window. Data reduction similar to that 
used in atmospheric neutrino analyses at SK3'0 is applied to select fully contained (FC) neutrino 
interactions. The criteria are: a) There is no detector activity within 30 /us before the event, b) 
The total collected p.e. in a 300 ns time window is > 200 («20MeV deposited energy), c) The 
number of PMTs in the largest hit cluster in the outer-detector is < 10. d) The deposited energy 
is > 30 MeV. Finally, a fiducial cut is applied accepting only events with fitted vertices inside 
the same 22.5kt volume used for SK atmospheric neutrino analysis. The detection efficiency of 
SK is 93% for CC interactions and 68% for NCi ne i interactions, for a total of 79%. Similarly to 
the lkt, the inefficiency is mainly due to the energy cut. Figure || shows the A(T) distribution at 
various stages of the reduction. A clear peak in time with the neutrino beam from the KEK-PS 
is observed in the analysis time window. 

In order to compare the expected and observed E v spectrum, the neutrino energy was 
reconstructed for the 1-ring //-like FC events in the SK fiducial volume assuming charged-current 
quasi-elastic interaction {y + n — * fi + p). For this interaction, the initial neutrino energy can 
be calculated from the momentum^) and angle(#) of the outgoing muons as 

m N E^ - ml/2 
m N - E^+p^ cos 9fj, ' 

where rriM and are masses of nucleon and muon, respectively. It should be noted that for non- 
quasi-elastic interaction events, the neutrino energy is underestimated. The expected spectrum 
was made from the beam and neutrino interaction MC simulation. Overall normalization was 
done using the event rate measurement at the lkt. 

3 Results 

Fifty-six FC events are observed in the fiducial volume. The 1.5/iS selection window gives an 
expected background from atmospheric neutrino interactions of order 10 -3 events. 



(1) 
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Figure 3: Far site timing analysis: (top) A(T) distribution over a ±500/is window for events remaining after cuts 
a) (unfilled histogram) and b) (hatched), and after all cuts except timing (shaded). The cuts are described in the 
text, (bottom) A(T) for a ±5/j,s window for FC events with vertices in the fiducial volume. 



The systematic uncertainty of the lkt measurement is 4%, for which the leading terms are due 
to vertex fitting and its effect on the fiducial cut. Other sources such as energy scale uncertainty 
and the treatment of multiple events in a spill are relatively small. The statistical uncertainty 
of the lkt measurement is < 1%. The systematic uncertainty of the extrapolation from the near 
site measurement is estimated to be i 7 % based on the PIMON measurement uncertainty and 
beam simulation uncertainty for low energy neutrinos. The systematic uncertainty in the SK 
measurement is 3%, mainly due to the fiducial cut. The systematic uncertainty term coming 
from uncertainty in the neutrino energy spectrum and cross section is small due to cancellations. 
The quadrature sum of all known uncertainties for the far site event rate prediction is 1liQ%- The 
resulting expectation is 80.6^ 8 'q events in the absence of neutrino oscillations. The observed and 
expected numbers for various categories are summarized in Table |]. The number of Cherenkoy 
rings and particle identification are reconstructed by the same algorithms as those used at Srfl, 
and the event category definitions are also the same. The expected number of events for analyses 
of CC interactions in the other near detectors are 8Tt]j[ for MRD (iron target) and 88^3, for 
SciFi (water target jH, each of which is consistent with the value based on lkt analysis. 

Table 1: Summary of the observed and expected numbers of FC events in the fiducial volume at the far site. 



Event Category 


Observed 


Expected 


Single Ring //-like 


30 


44.0±6.8 


Single Ring e-like 


2 


4.4±1.7 


Multi Ring 


24 


32.2±5.3 


TOTAL 


56 


80.6t^ 



Figure |] shows the reconstructed E u spectrum for the 1-ring //-like FC events in the SK 
fiducial volume. The histogram in the figure is the one expected from the beam and neutrino 
interaction MC. The difference between the expected and observed spectrum indicates the dis- 
tortion due to the neutrino oscillation. 

4 Conclusion and Future Prospect 

The crucial principles of K2K have been proven to work: The beam direction has been contin- 
uously monitored and controlled within ±lmr for several months duration by controlling the 
proton position on target with millimeter accuracy and operating the horn magnet stably. Pion 
kinematics have successfully been measured in situ allowing prediction of the neutrino beam at 
the distance of 250 km with six to seven percent accuracy. Events at the far site have been iden- 
tified by means of GPS timing, reducing backgrounds to a negligible level. Fifty-six FC neutrino 
events have been observed where 80.6Zts'o are expected. The study of the spectral shape for the 
oscillation analysis is now on going. In order to confirm the neutrino oscillation phenomena and 
to determine the oscillation parameter precisely, the beam spectral shape, neutrino interaction 
and their errors are being carefully studied. For more complete study of the beam spectrum and 
neutrino interaction, a new fine-grained (~2cm) detector will be installed in 2003. It consists 
of 3m long plastic scintillator bars and wavelength shifter fibers. The experiment expects to 
accumulate 10 20 protons on target, providing sufficient statistics to study neutrino oscillations 
by spectral analysis for disappearance. 
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Figure 4: Spectrum of reconstructed E u for the 1-ring /i-like FC events in the SK fiducial volume. The histogram 
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